Quantum bath refrigeration towards absolute zero: unattainability principle 

challenged. 
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A minimal model of a quantum refrigerator (QR), i.e. a periodically phase- flipped two- level 
system permanently coupled to a finite-capacity bath (cold bath) and an infinite heat dump (hot 
bath), is introduced and used to investigate the cooling of the cold bath towards the absolute zero 
(T — 0). Remarkably, the temperature scaling of the cold-bath cooling rate reveals that it does not 
vanish as T — > for certain realistic quantized baths, e.g. phonons in strongly disordered media 
(fractons) or quantized spin- waves in ferromagnets (magnons). This result challenges Nernst's third- 
law formulation known as the unattainability principle. 



Introduction. One of the generally unsettled funda- 
mental problems of thermodynamics is the nature of the 
ultimate limitations on cooling to absolute zero, T = 0. 
Attaining T = in a finite number of steps, or, more gen- 
erally, in finite time, is prohibited by Nernst's unattain- 
ability principle which is the dynamical formulation of 
the third law of thermodynamics [lJ-IH .However, the uni- 
versality of this principle has been postulated rather than 
proven. It is also debatable whether this formulation is 
always equivalent to Nernst's heat theorem, whereby the 
entropy vanishes at T = O.Do both formulations of the 
third law hold for all quantum scenarios? The investi- 
gation of this open fundamental problem, pertaining to 
quantum refrigerator (QR) schemes [4-8], raises several 
principal questions: (i) How does the cooling rate scale 
with the bath temperature and does it necessarily vanish 
as T — >• 0? (ii) Does a QR differ from its classical coun- 
terpart regarding compliance with the third law and to 
what extent is such compliance model-dependent? An- 
swers to these questions are important not only to the 
understanding of the foundations of quantum thermody- 
namics 0, but also to the design of novel QR schemes 
compatible with the needs of quantum nanotcchnologies 
that require compact (nanosize) coolers capable of ultra- 
fast cooling (lo| . 

Here, we propose and explore the simplest QR de- 
sign thus far that allows us to address the fundamen- 
tal issues raised above. The working medium is a single 
two-level system (qubit), permanently (rather than inter- 
mittently, as is done in traditional cycles coupled 
to a finite-capacity bath to be cooled and to another, 
much larger and hotter, heat dump. The pumping oper- 
ation consists of fast modulation of the qubit energy by 
means of periodic 7r-flips of the qubit phase. We find that 
this QR can cool down a finite-capacity (yet macroscopic 
and spectrally-continuous) bath only if the modulation- 
period is within the bath-memory (non-Markovian) time. 
Hence, the cold-bath spectrum is crucial in determining 
the cooling condition and rate. Our most striking finding 
is that for certain experimentally realizable baths, such 



as quantized spin-waves in ferromagnets (magnons) 
or acoustic phonons in strongly disordered media (frac- 
tons) [lH, the cooling rate remains finite as T — > 0, in 
apparent violation of the dynamical formulation of the 
third law. 

Model and analysis. A control qubit is weakly cou- 
pled to two baths via the system-bath interaction hamil- 
tionan: Hsb = ct x (Bh + Be), where o~ x is the spinor 
x-componentj-Bc 1 is the operator of a finite cold bath ( 
C) which we wish to refrigerate, and Bh that of a much 
larger hot bath ( H) that remains nearly unchanged. The 
qubit energy is periodically modulated by an external 
field v(t) via the Hamiltonian H ext = ^a z v{t). An il- 
lustration (Fig. 1-inset) is that of a charged quantum 
particle in a double-well potential that is periodically 
phase-flipped by off-resonant pulses and is coupled to a 
spatially-confined (macroscopic) C-bath to be cooled, as 
well as to a nearly-infinite H environment into which the 
heat is dumped. This scheme bears analo gy t o radiative 
(sideband) cooling in solids and molecules [13j, if one vi- 
sualizes the red- and blue- shifted qubit frequencies as 
Stokes and anti-Stokes lines respectively. 

Our analysis reveals the crucial role of the quantized 
characteristics of system-bath coupling in determining 
the attainability of T — > 0. By contrast, the results are 
insensitive to the QR scheme chosen (see Discussion). 

The general condition for steady-state refrigeration, 
under periodic, off resonant, modulation, is positive heat 
current from C to H via the qubit. The sign and mag- 
nitude of the current is determined by the steady-state 
solution of a non-Markovian master equation (ME) for 
the qubit density matrix [14J. The ME, which is ac- 
curate to second order in the system-bath coupling, al- 
lows for time-dependent modulation of the system that 
is much faster than the bath-memory time t c . It is valid 
at any T, as shown both theoretically [l4[ and experi- 
mentally [15l.ll6 |. Deviations of the evolution (Suppl. A) 
from that described by the non-Markovian ME include 
system-bath entanglement (correlations) effects that can 
be compensated by readjusting the qubit excitation, as 
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well as bath dynamics effects (violation of the Born ap- 
proximation whereby the bath is constantly in a thermal 
state) [l7|]- Yet such deviations are of fourth-order in the 
system-bath coupling and thus negligible for weak cou- 
pling. The cooling of a finite-capacity bath is the result 
of infinitesimal temperature changes over many modu- 
lation cycles, consistently with the Born approximation 
underlying the ME. Since the Born app roximation is the 
more accurate the larger the bath[18|, we assume that 
the finite-capacity bath is macroscopic and has a contin- 
uous spectrum, which does not exhibit mode discreteness 
or recurrences that may otherwise invalidate this approx- 
imation and bath thermalization altogether flj| . 

Only the diagonal elements of the qubit's density ma- 
trix ps (energy-state populations) play a part here, al- 
though the ME also allows for coherences (off-diagonal 
elements) [14], but these are absent at t = (starting at 
equilibrium) and remain so under the modulation. The 
quantumness of the ME, even when it is diagonal in the 
energy basis, is embodied by the qubit interlevel tran- 
sition rate and their non-Markovian time-dependence 
(Suppl. B). Periodic phase shifts of the qubit at intervals 
r dynamically control its coupling to the baths and the 
resulting transition rates. When r is comparable to the 
bath memory-time t c , these phase shifts modify the de- 
tailed balance of the transition rates and thereby either 
heat or cool the qubit depending on r [HI, E3] • In what 
follows, we analyze the steady state and the slow changes 
of the bath temperature as a result of these periodic per- 
turbations. 

Under weak-coupling conditions, the qubit evolution 
caused by the baths is much slower than r ~ t c . Hence, in 
steady state, we can use time-averaged level populations 
and transition rates between the periodically-perturbed 
qubit levels (Suppl. B). These time-averaged (steady- 
state) equations can be recast, upon introducing the po- 
larization of the qubit S = (p ee — p gg )/2, into 
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(1) 



Here the |e) — > \g) averaged transition rate from the ex- 
cited (e) to the ground (g) state, is R e and its \g) — > |e) 
counterpart is R g . The averaged transition rates for 
t ^> r are found, upon expanding the qubit energy un- 
der periodic frequency modulation v(t) into the harmonic 
(Floquet) series (Suppl. B, [Igj]) 



R e(g) =2irJ2 P ™ G T[Huo+mA)]; 
m 

Prn = \e m \ 2 , S m = ~ f e* X ("(*')-"o)dt' e imAt dt, (2) 



Here m are all (positive and negative) integers, P m arc 
the probabilities of shifting Gr(w) by mA, A = — , 



from the qubit average frequency aj , Gt{oj) being the 
temperature-dependent bath-coupling spectrum, i.e. the 
Fourier transform of the bath autocorrelation function: 
G t {uj) = e luJt (B(t)B(0))dt = e"/ kBT G T {-u). For 
a bosonic bath: 
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GtH =Gf (w)(n(w) + 1); G (w) = \ 9 (uj)\ 2 p(uj); 
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(3) 



g(us) being the system-bath coupling, p(uj) the bath - 
mode density and n(u>) the w-mode thermal occupancy. 
These expressions are also obtainable by Floquet (har- 
monic) expansion of the periodically-driven Markovian 
(Lindblad) Liouvillian [H,^. 

In the presence of hot (H) and cold (C) baths, under 
the assumption of a qubit weakly coupled to both baths, 
the transition rates are split into additive harmonic con- 
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Eq.(l) is also split into (Suppl. B) S = J2 m ( S m + S m), 

where S^ H ^ is the m-harmonic polarization flow caused 
by the cold (hot) bath only. The averaged heat flow,Q, 
through the qubit is correspondingly divided into the C 
and H bath-contributions. The steady-state Eq. (1) then 
gives rise to 



Jc{H) = Qc(H) =^2(^0 + mA)S. 



C(H) 



(4) 



which is the sum of rates of heat-exchange with the 
respective baths via wo + ?nA quanta. Positive Jc implies 
refrigeration, i.e., heat flow from the cold bath to the hot 
bath via the modulated qubit . 

It is advantageous to use periodic, alternating, 7r-phasc 
shifts (phase flips) as they give rise, to leading order, 
to two symmetrically opposite frequency shifts of Gt 
at wo ± A. Then, by the Floquet expansion we ob- 
tain the probability distribution wherein Pq = and 
(2/7r) 2 are the leading terms 141.121 [. 
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Let us choose sufficiently large A, of the order of the 
spectral width T = l/t c , which is the inverse memory 
time of the cold bath, such that at uj ~ w + A the qubit 
is coupled only to the hot bath, while at uj ~ loq — A 
it is coupled to both the cold and the hot baths. More 
precisely, we require that 

Gr(w + A) w G%(lo + A) > G%(lo + A); 
G% (uj + A) > G$ (cj - A) , G% (w„ - A) (5) 

where G^S H \lo) is the respective temperature- 
dependent bath-coupling spectrum. This requirement 
can be satisfied if the cold bath (C) is spectrally localized 
with upper cutoff u> cu t < ujq + A. By contrast, for the hot 
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FIG. 1. Main panel: Schematic depiction of the required bath 
spectra and the qubit frequency shifts due to periodic phase 
flips. Inset: Schematic realization of the modulated qubit and 
its coupling to the baths. 



pump requires the qubit to be simultaneously coupled to 
the C and H baths at loq — A and ujq + A, respectively. 

Cooling rate scaling with temperature. In what follows, 
we investigate the QR action (heat pumping from C to 
H) under the assumptions that the hot bath is practically 
infinite, hence Tjj = const, whereas the macroscopic cold 
bath has finite heat capacity, cy < oc, resulting in slow 
evolution of Tc it) under the QR action. To estimate this 
evolution we use the standard thermodynamic definition 

M 



dT c (t) 

Cv ^T = Jc = 



(8) 



which presumes that Tc is well-defined at all t (since the 
bath has a continuous spectrum and is large enough to 
thcrmalizc at finite times). 

In order to infer the temperature-dependence of the 
cooling rate we shall examine the scaling of cy and 
Jc with Tc : a) The constant-volume heat capacity of 
the cold bath, cy, depends on the dimensionality of the 
bosonic bath. If p(u>) ~ oj d ~ 1 is the d-dimensional density 
of modes and Tc <C io cu t (fcs = h = 1), then 



bath (H), the required rise of Gtf with oj is obtained for 
most common bath spectra, provided the cutoff of (w) 
is much higher than u cut of G^(w): e.g., for phonons in 
bulk media or photons in open space, (u>) oc oj 3 satis- 
fies Eq. ©. 

Under the conditions of Eq. (JSJ) we find that the 
steady-state heat current from C to H is (Suppl. B) 



J c = (wo ~ A)S% = 



(wo-A)- 



(6) 

A)[n c (w - A) - n H {uj + A)] 



[2n H (uJo + A) + l] 



The balance of Jc and Jw (cold and hot ) currents fl3J) 
obeys the second law 0, [22[: It can be verified that the 

entropy production rate ^ satisfies: ^"(f^ + j^j > 
for any initial state. 

From Eq. (6), the heat pump (QR) condition Jc > 
amounts to 



n c (uj - A) > n 11 (uj + A) <*=> 



An analogous relation holds if n c ^ H \ui) are Boltzmann 
rather than Bose factors (occupancies). 

Equation (|7|) reveals the crux of the heat pumping 
(QR) effect: although by definition n c (uj n ) < n H (coo), 
heat can flow from the cold to the hot bath if the C-bath 
thermal occupancy at luq — A is higher than that of the H- 
bath at u ~ ujq + A. If A is too small for Eq. ([7|) to hold, 
we recover the natural heat flow direction H — >• S — >• C 
at steady-state. In addition, Eq.© implies that the heat 
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(9) 



b) The scaling of the cold-bath heat current, Jc, in 
Eq. ||7J) can be deduced if we maximize the heat flow [24[ 
with respect to A (our control parameter). This gives 
the dependence of ujq — A ~ Tc @, Hence, to main- 
tain the maximum heat flow, we have to slowly increase 
A with time, so as to approach Tc —> 0. The closer 
to Tc —> 0, the lower is ljo — A, hence the steady-state 
dynamics |(7J) and its slow change ([5J become increas- 
ingly more accurate. Correspondingly, we parametrize 
Jc in Eq. (6) using Eq. (3) and assuming the low- 
frequency range of the cold bath < uj = ujq — A u cu t'. 
lim^o |ff(u;)| 2 oc w 7 , p(u>) w u^ 1 . Here, |.g(a;)| 2 is 
the 7- dependent system-coupling to the bosonic bath 
(discussed below). The heat current, maximized for 
ujq — A w Tc, then obeys the scaling 



J C (T C ) cx 



(10) 



c) Upon substituting Eqs. l[9])-([T0]l in Eq. © we ob- 
serve that the T c scaling of cy is canceled by a similar 
scaling of the density of modes in Eq. (JTUJ) . The resulting 
scaling yields 



dTc/dt 



-ATI. 



(11) 



Here the constant A oc 1/V: the larger the bath the 
slower its cooling. 



4 



Remarkably, the scaling only depends on the 7- 
th scaling power of the system-bath coupling strength 
|g(w)| 2 . For different 7 the time-dependence of Tc, 
starting from the same Tc (0), is plotted in Fig. [2] 
For 7=1 we have exponentially slow convergence to 
Tc — > 0, conforming to the third law. Yet, strikingly, for 
< 7 < 1, Tc(t) — > at finite time, thus violating the 
accepted dynamical formulation of the third law [l|-|3j , if 
the frequency-dependent coupling |g(w)| 2 is sub-linear hi 

In what follows, wc examine the possibility of such 
scaling for different bosonic baths. To this end, consider 
a qubit immersed in a periodic medium, whose local dis- 
placement is a linear combination of normal-mode cre- 
ation and annihilation operators (bath excitations/de- 
excitations) B(x) = ^X; c -^=(^ c (x)at(k) + h.c). 

The normal-mode functions are labeled by the wave 
vectors k that belong to a reciprocal lattice bounded 
by the Debye cutoff (w(k) < Lu cut = Wd). The cou- 
plings of a charged or dipolar system to bath excita- 
tions/deexcitations are to leading order determined by 
the gradient of the displacement operator 

VS(x) = ^^^l=(V^(x) a t(k)_/ l . c .) (12) 

vv c \Mk) 



When 0g(x) = 
stant scales as 



e lk x the corresponding coupling con- 
<?(w(k))| ~ k m . We can discern three 



generic types of scaling of the coupling constant: 

i) For acoustic phonons w(k) ~ v\k\, where v is a sound 
velocity and the coupling strength satisfies |g(w)| 2 ~ w, 
i.e. 7=1. Therefore, acoustic phonons used as a cold 
bath do not violate the dynamical third law formulation: 
the optimal cooling to zero temperature is exponential in 
time. 

ii) Amorphous (glass) materials may exhibit effects of 
fractal disorder. These effects imply different scaling of 
the displacement of the mode function 0g(x), |V</>g(x)| ~ 
w 7 |0g(x)|: normal phonons are replaced by fractons for 
which 7 takes fractional values. In particular, for some 
materials 7 < 1 [l^ . Hence, for a cold bath composed of 
such fractons the violation of the third law is expected. 

iii) Another system which leads to a violation of the 
third law is the magnon (spin-wave) bath in a fer- 
romagnetic spin lattice with nearest neighbor interac- 
tions, below the critical temperature. The Holstcin- 
Primakoff transformation of the jth spin Pauli matrix 

g+ _ a , .-a /OOM/2/1 _t_ /no\V2. 
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(2S)^(l-a] aj /2Sy 



to boson annihilation and creation operators a 



al- 



lows to represent the system as a set of interacting 
harmonic oscillators. Introducing the collective spin- 
wave (magnon) variables a(k),er(k) satisfying aj = 

^= e _lk x 3a(k),wc can rewrite its Hamiltonian in 



the form Hq — ^gu;(k)a^(k)a(k) + Higher order terms. 
At low temperatures the nonlincarity in the Holstcin- 
Primakoff transformation can be neglected and the sys- 
tem becomes equivalent to a bosonic system governed 
by the Hamiltonian Hq, whereby the dispersion law is 
quadratic in the low- frequency region, w(k) ~ (|k| 2 + 
constant) . The local spin variable aj can then be directly 
coupled to the qubit by a dipole-dipole (spin-spin) inter- 
action. Hence, the main difference between the dipolar 
coupling to acoustic phonons and magnons is the absence 
of the dispersive-coupling coefficient , k _ for the lat- 

ter. Therefore, the coupling strength to magnons satis- 
fies |g(w)| 2 ~ 1 (7 = 0), which implies the violation of 
the third law for magnons. 

Discussion. We have analyzed the cooling process of 
a bosonic bath towards the absolute zero using a new 
minimal model of a quantum refrigerator: a single two- 
level system (qubit) permanently coupled to a spectrally- 
restricted cold bath with finite heat capacity and a hot 
bath with infinite heat capacity has been shown to act as 
a heat pump, under appropriate modulation. The heat 
flow is proportional to the population-difference of a pair 
of oppositely-shifted bath modes that are selected by the 
qubit modulation (phase-flip) rate, analogously to side- 
band cooling [I3J . The attainable cooling rate challenges 
the third law of thermodynamics, in the sense that arbi- 
trarily low temperature of the cold bath may be reached 
in finite time by the heat pump for certain quantized cold- 
bath spectra: e.g. magnon and fracton baths. 

In solid-state ferromagncts or glasses, interactions of 
control qubits with other baths unaccounted by the 
model, as well as tiny deviations from the predicted weak- 
coupling, steady-state dynamics (discussed in Suppl. A,B 
) may restore the third law. Nevertheless, surprisingly 
fast cooling (7 < 1) may still be observed down to some 
(material-dependent) temperature. It would be prefer- 
able to demonstrate this effect for quantum dots cou- 
pled to controllable baths composed of nuclear spins in 
solids [2^ or for atomic dipoles in optical lattices 27 1: 
in both cases the systems are highly shielded from other 
baths, while the lattices can be engineered to conform to 
the nearest-neighbor ferromagnetic model that engenders 
magnons. 

This study of the compliance with the third law for 
quantized non-Markovian baths indicates that the tem- 
perature scaling of the cooling rate is not specific to the 
chosen QR model; it is similar to the scaling obtained 
for the very different noise -driven QR[25|]. Namely, the 
scaling is not sensitive to the form of driving, nor to the 
method of treating the steady-state dynamics. Hence, 
the dependence of the scaling on the system-bath cou- 
pling dispersion is general. It provides new insights into 
the bounds of bath cooling in quantum thermodynam- 
ics. It shows that Nernst's principle of unattainability 
of the absolute zero in finite time may fail and is 
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1 2 3 4 5 6 

FIG. 2. Tc change with time (cooling) for three different 
system-bath coupling-strength dispersion laws: 7 = 1 (acous- 
tic phonons), 7 = 3/4 (fractons), 7 = (magnons). 

not always equivalent to Nernst's heat theorem (see In- 
troduction): the latter holds true since a bosonic bath 
has a unique ground state whose entropy must vanish. 
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SUPPLEMENTARY A: UNIVERSALITY AND 
ACCURACY OF THE NON-MARKOVIAN 
MASTER EQUATION 

The dynamics derived in the main text can be for- 
mulated for the following Hamiltonian for a phase- 
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modulated qubit (TLS) weakly coupled to a general bath: 

(Al) 
(A2) 

H SB = a x B, (A3) 



H t ot = Hs + Hb + eHsB , 



where Hb is any time-independent bath Hamiltonian and 
B is any time- independent bath operator and 



max(?7 fe /w ), 



(A4) 



is a small dimensionless parameter normalizing the rate 
rjk of the maximally coupled bath mode divided by the 
TLS resonant frequency, u>o- 

As shown in Ref.[3l, one can derive, to second order in 
the system-bath coupling, the following non-Markovian 
master-equation for the reduced system density matrix, 
Ps(t): 



Ps(t) 



i[H s ,ps(t)} + 
S(t,T)p s (t),a x 



H.c 



where 



e 2 (Be- lHBt Be iHBt ) 



$ T (f) 

is the correlation function of the bath and 
S(t,r) 



e -iHs(t-T) a e iH s (t-r) 



•} (A5) 



(A6) 



(A7) 



If we start at equilibrium, where ps is diagonal in the 
energy basis of the TLS (|e), \g}) then it remains so un- 
der the action of the diagonal Hs(T). This ME does 
not invoke the rotation- wave approximation [14 , Il7| and 
hence allows for arbitrarily fast modulations of the sys- 
tem, which cause its anomalous heating or cooling on 
time scales comparable to lo^ 1 . The resulting rate equa- 
tions are then given by 



Pee(t) = -f>gg(t) = R g (t)p gg - R e {t)p e 



(A8) 



where the non-Markovian time-dependent rates R g { e ) (t) 
are given for periodic phase flips in Suppl. B. It yields 
a Floquet expansion of p ee (t) and the corresponding po- 
larization, as detailed in Suppl. B. 



Accuracy of the Non-Markovian Master Equation 

The second-order (Zwanzig-Nakajima type [Hj]) non- 
Markov Master Equation (ME) involves two crucial ap- 
proximations: The first is that the initial state is a prod- 
uct state of the qubit and bath density matrices. The 



second is that the bath state ps does not change during 
the evolution (Born approximation). These approxima- 
tions yield: 



Ptot{t) 



ME 



Ps(t) <8) PB 



ME 



(A9) 



There are two types of deviations of the ME-described 
initial state from its exact (initial state) evolution. 

1. The first results from the system- bath correlation 
effect that does not enter the second-order ME cal- 
culation. This deviation raises the effective qubit 
temperature above that of the bath, by an extra 
excitation, 

S Pee = \{{a z {e)) - (a z (e = 0))) ~ 0(e 2 ). (A10) 

It can be compensated upon using the exact state 
populations [l7j : 



Pee(t) = Pee{t) + 5 Pee, Pgg{t) = 1 - p ee (i)(All) 

ME 



This correction vanishes over time intervals > 
and does not affect the steady-state polarization 
(heat-flow) rates considered in Suppl. B. 

2. The other type of deviations is due to the bath: the 
ME assumes that ps is in a Gibbs (thermal) state, 
while the exact state has small deviations from it, 
responsible for correlations among the bath modes. 
These are, in general, complicated functional forms 
of temperature and frequencies (see Suppl. Ma- 
terial of [It}, particularly Gordon et al., (NJP 11, 
123025(2009)). The deviation from the Gibbs state 
used in (|A9[) is only of fourth order in e, which can 
be neglected in the weak-coupling limit. 

More generally, we can show that the deviations from 
master equation analysis are always to fourth order in 
e. To this end, the total density matrix (of the sys- 
tem+bath) at any time t can be expanded in powers of 
e as 



p tot (t)=pi°l(t) + epill(t) + eipill(t) 



■ (A12) 



where p\l t (t) need not necessarily be a valid density ma- 
trix. Upon inserting (|A12[) in the equation of motion for 
Ptot(t), and collecting terms of successive orders of e, we 
get 



Ps 



ieTr B [H SB {t),pto 



-I dt'Tr B [Hs B {t),[H SB {t'),p tot {t)}}+0{e% (A13) 
'o 
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where the time-dependence of is determined by 

Hs(t). It is important to note that odd powers of e do not 
contribute, since the trace of the interaction Hamiltonian 
Hsb vanishes. 



One can see upon collecting the even powers of e in 
Eqs. (|Al2j) . (|Al3j) . that the corrections to the master- 
equation analysis are always of fourth order in e. 



SUPPLEMENTARY B: FLOQUET-EXPANSION 
OF POLARIZATION AND TRANSITION RATES 



R e (t) = 2 Re / At' exp[iuj {t - t')]e(t)e*(t')<f> T (t - t') 



E< 

k,l 



k£i 



AluGt(ijj) 



+ f cos[(ojfc — u)i)t] sin[(wo + ~~ w )£] 



2 sin[(ojfc — sin 2 [(o;o + u>i — ui)t/2] 

LUq + tjJl — u> 

R g (t) = 2 Re f At' cxp[-io; (t - t')}e*{t)e(t')^ T (t - t') 
Jo 



= 2 



AujGt(ijj) 



E, f cos[(w fc - uji)t] sin[(w + u>i + uj)t] 



k-.i 



LJq + UJl + UJ 



2 sin[(ojfc — 0Ji)t] sin 2 [(a;o + u>i + uf)t/2] 



Under energy modulation and in the presence of weak 
coupling to two baths, the non-Markovian ME leads to 
the following additive contributions to the expressions for 
qubit polarization 



s = s° +s 

(rC + rC)S+^- 



t>C tdH _ r>H 

-^-« + i?f)S+ g o e . 



(Bl) 



Under impulsive 7r-flips of the qubit phase at periods 
t, the Floquet expansion of e(t) yields £o = 0, £±1 = 
±2/tt, e k> i = and uj k = (2k + l)n/r. These are 

respectively the fcth-harmonic amplitude and frequency 
of the periodic 7r-flip modulations. 

Taking into account the additive contributions of cold 
and hot bats to G T (u), i.e. G t (uj) = G$(uj) + G"(uj) 
and assuming the conditions of Eq.(5) of the main text, 
we obtain the time-averaged partial decay rates as 



R, 



G^u:o-A)[n c (u Q -A) + l}, 



From Ref. [14| we have the |e) — > \g) and \g) |e) 
transition rates, respectively, as 

R e (t) = 2 Re / At' exp[iw (i - t')}e(t)e*(t')^ T (t - t'), 
Jo 



R g (t) = 2 Re/ At'exp[-iuo(t-t')}s*(t)E(t')<f> T (t-t'), 
Jo 

(B2) 

Here "Re" stands for real part, e(t) for the modulated 
phase factor (a unimodular periodic complex function) 
ujo is the qubit's resonance frequency and the bath- 
correlation function = / AgjGt{w) cxp(— iu>t) is 
the Fourier transform of to the bath coupling spectrum 
G T (u). 



The Floquet expansion of e(t) can be substituted into 
Eq. (|B2|) . to obtain (k, I being integers) 



R 

R, 
R 



H 



-_Gq(^o ~ A)n c (ujo - A), 
[G$(u; + A){n H (u; + A) + l] 
-G^K + A)n H ( Wo + A). 



(B3) 



where A = 2-k/t. Under these assumptions, we can 
find a periodic quasi-steady state with constant polariza- 
tion Ss , namely, with the use of Eq. (|B3I) 



S s = - 

G%(u -A)+Gg (uq+Au) 



Using this result and separating the C and H contribu- 
tions, we obtain the heat flows in the steady state 

G^C^o-AjG^ (ajQ+A)[ng(qjQ-A)-ng(^o+A)] 



G<f (uo-A)[2n<£(u„-A) + l]+Gg (w +A)[2n»(w +A) + l] 

- 



(B5) 
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For arbitrarily wide C and H coupling spectra, the heat 
flow reads 



AC(m) 







(w + 


A)[n c {uj 


-A) 


- {ojq 


+ A)] 




+ A)G " 


(w - 


K 

A)[n c (oj 


+ A) 


- n§ (uj 


-A)] 




+ A)G " 


(wo + 


K 

A)[n c (w 


+ A) 


- (loo 


+ A)] 


G{f(wo 


-A)G ff 


(w - 


K 

A)[n c (w 


-A) 


- n% (u> 


-A)] 



K 



(B6) 



where 



K = G$(uj - A)[2n£(w - A) + 1] 
+ Gq 7 (wo + A)[2n£(w + A) + 1] 
+ G"(ujq - A)[2n%(io - A) + 1] 
+ Gq(ojq + A) [2n£ (w + A) + 1] . (B7) 

In the case of flat and wide (Markovian-like) spectra, all 

four terms are of comparable magnitude, hence Sg is 
negative, i.e. there is no QR action. By contrast, under 
the assumption of Eq.(5) of the main test, the last three 
terms are negligible, yielding Eq.(7) and possible change 
of sign of the heat flow for QR. 



Steady states and currents by Floquet expansion of 
the Lindblad operator 

An alternative method for performing the Floquet ex- 
pansion of the steady-state of a periodiaclly- flipped qubit 



coupled to two baths is based on the Lindblad operator, 
as expounded in the tutorial in Ref. fl9| 

For the qubit density p to be a steady state of the 
Lindblad operator we expand C — ^2C 3 where j = H, C, 
and q is the Floquet harmonic. The expansion yields (see 
notation in the main text) 



C ip=^lf (°J N + ?A) ( [a- p, a+] + [a- , pa+] ) + 
Gj(-w - qA)([*+p,*-] + [a+,pa-])) (B8) 

The qubit steady-state then has the form 

A 
B 



A E qj P(q)G 1 (Lo + qA)e 



"0 + <l& 



(B9) 



by 



where A = — . The cold (hot) current is then given 



J c(H) = ^?(9)K+?A)G c < ff )( W o+9A)! 



'O(H) 



A 
B 



- 1 

(BIO) 



The magnitudes and signs of these steady-state cur- 
rents are the same as for the time-averaged ME solutions 
detailed above. 



